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PREFACE 

This  work  on  the  equation  of  state  was  done  initially  in 
1967.  It  has  led  to  an  equation  of  state  'or  pure  water  based  on 
the  work  of  Kell  and  Whalley.  Subsequent  work  by  Millero.t  in- 
cluding corrections  in  the  data  of  Kell  and  Whalley  based  on  his 
results  have  yielded  a more  complex  equation  of  state  for  pure 
water  and  sea  water.  Tlie  results  in  this  report  are  useful  not 
only  in  developing  an  easily  differentiated  equation  of  state  of  sea 
water  suitable  for  a number  of  purposes  requiring  less  accuracy 
than  Millero’s  equation  but  providing  a basis  for  the  represen- 
tation of  the  PVT  properties  of  other  solvents  and  solutions  con- 
taining water  such  as  for  dioxane -water  mixtures. 

F.H,  Fisher,  October  16,  1975 

tF.J.  Millero,  R.A,  Fine  and  D.P,  Wang  "The  Squatlm  of  State 
of  Seawater, " J.  Mar.  Res.,  32,  pp.  433-456,  (1974). 
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ABSTRACT 


For  representing  the  PVT  properties  of  pure  water  and  solutions,  sea 
water  in  particular,  we  have  selected  the  Tumlirz  equation 


V = V®  - KjS  + 


X 

P0  + K2S  + P. 


where  V is  the  specific  volume  as  a function  of  pressure  P,  Salinity  S and 
temperature.  For  pure  water  S + 0 and  for  sea  water  V®,  X.  and  P0  are 
the  pure  water  values.  For  pure  water  we  have  obtained  an  analytic  func- 
tion which  tits  the  data  of  Kell  and  Whalley  to~8  ppm  from  0 - 100°  and- 
15  ppm  from  100  - 150°  up  to  1000  bars.  The  calculated  density  maximum 
is  at  4.00°.  Tlte  form  of  the  equation  appears  to  be  useful  not  only  for  rep- 
resenting the  properties  of  sea  water  but  for  deriving  partial  molal  prop- 
erties of  solutions  at  low  concentrations  from  high  concentration  data. 


INTRODUCTION 

The  goal  of  an  equation  of  state  of  a liquid 
is  to  represent  PVT  data  in  a simple  form  and  pro- 
vide insight  into  the  physical  properties  of  the  liquid. 
In  particular,  an  equation  of  state  for  a solution  such 
as  seawater  should  reduce  totliat  of  die  pure  solvent 
as  the  solute  concentration  goes  to  zero.  In  order 
to  understand  the  effects  of  adding  solute,  the  equa- 
tion of  state  for  the  solvent  should  be  established  as 
accurately  us  possible.  In  this  paper  we  consider 
two  equations  which  have  been  used  before  to  rep- 
resent electrolyte  solution  and  sea  water  as  well  as 
pure  water. 


Since  Eckart’s  ^/review  on  the  equation  of 
state  for  pure  water  and  sea  water,  further  mea- 
surements have  been  made  on  both  pure  water z/and 
sea  water  3>4»5/.  The  high  precision  specific  vol- 
ume data  Kell  and  Whalley  obtained  for  pure  water 
as  a function  of  temperature  and  pressure  have  been 
particularly  useful. 

Two  equations  have  been  used  for  fitting 
specific  volume  data  as  a function  of  pressure,  the 
original  Tait  equation^/ 


AV  A 

VjAP  TT+P 


(1) 
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and  the  usual  Tait  equation  (UTE)  as  McDonak!-/ 
denotes  it  in  his  recent  review  on  equations  ofsta’e, 

vP=vl  + Clog|±|  ,2) 

The  original  Tait  equation  is  equivalent  to  the 
Tbmlirz  equation- 


that  Eckart  used  in  his  analysis.  Non-linear  fitting 
techniques  for  this  equation  were  developed(see Ap- 
pendix I)  to  avoid  eirors  inherent  in  less  direct 
methods. 

We  find  that  Keil  and  Whalley  data  can  be 
described  to  their  stated  precision  by  using  the 
Thmiirz  equation  at  low  temperatures  (0-90°)  and 
the  UTE  at  high  temperatures  (50-150°).  Older  pure 
water  data  were  not  sufficiently  accurate  to  select 
one  equation  in  nreference  to  the  other.  Since  the 
Tumlirz  equation  describes  the  Kell  and  Whalley 
data  within  us  stated  accuracy  over  the  entire  tern  - 
perature  range,  we  have  chosen  it  as  die  basis  for 
representing  P-V-Tproperties  of: -ire  water  and  sea 
water.  The  V®,  \ , ->nd  P0  parameters  can  be  des- 
cribed by  simple  polynomials  of  temperature,  com  - 
pleting  the  equation  of  state  for  pure  Water. 

The  Tumlirz  equation  was  fit  to  isotherms 
of  Wilson  and  Bradley  an  ! Newton  and  Kennedy  sea- 
water data  at  constant  sa.inity.  While  the  equation 
described  the  data  tc  its  precision,  the  parameters 
of  the  equation,  ' O*  k t and  P0  varied  erratically 
with  salinity  and  temperature.  We  found  that  by  con- 
straining X tc  the  same  value  as  for  pure  water  at 
that  temperature,  V®  and  PQ  become  linear  functions 
of  salinity.  This  procedure  yields  an  equation  to 
describe  the  P-V-T-S  properties  of  seawater,  which 
differs  in  form  from  that  used  by  Eettart  and  Wilson 
and  Bradley. 

CALCULATIONS 
A.  Pure  Water 

For  their  measurements,  Kell  and  Whalley 
(20  temperatures  up  to  150°  and  26  pressures  from 
5 to  1000  bars)  claim  accuracies  which  vary  from 
5 ppm  at  low  temperatures  ;0-50°)  and  pressures  to 
40  ppm  at  high  temperatu.  es  and  pressures.  Al- 
though Kell  and  Whalley-  published  both  their  ex- 
perimental and  smoothed  data,  we  use  their  experi- 
mental data  to  a void  any  systematic  errors  introduced 
by  theli  smoothing  procedure. 

Acomputer  routine  was  written  to  fit  equa- 
tions (2)  and  (3)  to  an  isotherm  of  water  data.  The 
routine  calculated  parameter  values  so  as  to  mini- 
mize the  quantity, 

, 1/2 

„ « [I  <VV/n]  (4) 


where  Vp'  is  die  volume  predicted  by  the  given  equa- 
tion at  pressure  P and  the  indicated  summation  is  at 
each  pressure  along  a single  isotherm;  this  quantity 
is  the  standard  deviation  of  observed  data  points  a- 
bont  tiie  fitted  curve,  sometimes  called  the  standard 
error  of  tstima.e.  Hie  details  of  how  these  calcu- 
lations were  carried  out  may  be  found  in  Appendix  1. 
Tables  1,  II,  and  HI  show  the  results  of  calculations 
fitting  these  equations  to  the  data  of  Amagat?/  and 
Kell  and  Whalley:?/  and  the  much  higher  pressure 
data  of  Vedam  and  Holton!?/.  Figure  1 shows  the 
temperature  dependence  of  the  various  parameters 
for  the  Tumlirz  equation  as  determined  from  the 
Kell  and  Whalley  data.  Figure  2 shows  the  temper- 
ature dependence  of  the  parameters  for  the  UTE  also 
fit  to  the  Kell  and  Whalley  data. 


TABLE  1 

Fits  of  Tumlirz  and  Usual  Tait 
Equations  to  Amagat  Data. 


Tumlirz  Equation 

Temp. , 
iegrees 

P0  (bars) 

Nbar  cc/ton 
x 10‘ ^ 

V,  cc/gm 

"g 

-e°o 

o -< 
u X 
b 

0 

5812.6 

t.7228 

0.70379 

11 

5 

6421.  1 

2.0164 

0.686C 

27 

10 

6228  4 

1.8584 

0.70193 

15 

15 

6322. 9 

1. 8742 

0.70449 

19 

20 

6569.  9 

1.9833 

0. 69989 

23 

30 

6831.0 

2. 1026 

0. 69741 

30 

40 

6888.2 

2. 1027 

0. 70244 

27 

50 

6787.  1 

2.0447 

0.71069 

25 

60 

6704.  8 

2.0126 

0.71670 

39 

70 

6189.5 

1,7609 

0.73802 

51 

80 

5493. 9 

1.4483 

0.76529 

19 

90 

5467.4 

1.4262 

0.77006 

21 

?CX) 

5357.9 

1.4694 

0.76901 

28 

198 

3328. 4 

1,0941 

0. 82994 

91 

Usual  Tait  Equation 

B (bars)  C.  bar  cc/gm 

Vj  cc/gn. 

tr  cc/gm 

x 106 

0 

2697.7 

0.  3180 

i.OOOH 

13 

5 

3002. L 

0.3392 

1.00000 

29 

10 

2893. 1 

0.3204 

1.00026 

15 

15 

2930.6 

0.3176 

1.00086 

22 

20 

3052.3 

0.3242 

1.00174 

24 

30 

3185.1 

0. 3298 

1.00428 

27 

40 

3205.5 

0.3283 

1.00767 

26 

50 

3155, 1 

0,3237 

1.01192 

23 

60 

3109.8 

0.3219 

1.01684 

34 

>n 

2850.8 

0.3032 

1.02249 

44 

8. 

2509.5 

0. 2789 

1.02889 

14 

90 

7453. 1 

0.2812 

1.  03574 

27 

100 

2449,8 

0.2904 

1.04322 

36 

198 

1425.3 

0.3298 

1. 15868 

82 

-2- 
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TABLE  11 


Fit  of  Tumlirz  anil  Usual  I’Jit  Equations  to  Kell  and  VVhalley  Data. 


Kell  and 
Whalley 
0.000 
10.002 

11.997 

24. 998 

24.998 
25.005 
30. 112 

39. 999 

50.007 
60.001 

70.003 

80.003 

90.007 


110.003 


1150.016 


Tumlirz  Equation 

Usual  Tait  Equation 

Pq  (bars)  \ bur  cc/gm 
x 10'1 

V,  cc/gm 

<r  x 10h 
cc/gm 

B bars 

C bars  cc/gm 

Vxcc/gm  3 

x 106cc/gm 

5825. 80 

1.73619 

0.702201 

4. 5 

2685.04 

0.31773 

1.000101 

11.9 

6323. 46 

1.92008 

0.696711 

6.2 

2932.58 

0.32553 

1.000315 

11.9 

6644.  11 

2.03391 

0.695724 

4.6 

3091.72 

0.32919 

1.001808 

9.0 

6722.73 

2.05508 

0. 697320 

6.1 

3131. 12 

0.32899 

1.002973 

10.3 

6723.00 

2.05585 

0. 697262 

5. 1 

3131.51 

0.32904 

1.002971 

9.0 

6729. 02 

2.05889 

0.697041 

5. 1 

3134. 15 

0.32931 

1.002974 

9.5 

6773. 51 

2.06654 

0. 699357 

3.5 

3155.79 

0. 32845 

1.004410 

6.5 

6783. 43 

2.05450 

0.705013 

4.3 

3160. 30 

0. 32605 

1.007847 

5.3 

6772. 43 

2.04980 

0. 709488 

4.9 

3155.27 

0.32584 

1.012119 

7.8 

6591. 20 

1. 96582 

0.718884 

5.0 

3064  46 

0. 32048 

1.017097 

5.2 

6433. 54 

1.90833 

0. 726164 

4.8 

2986.  09 

0.31825 

1.022749 

5.3 

6189, 63 

1.81601 

0. 735672 

6.0 

2864.  »1 

0. 31395 

1.029029 

4.1 

! 5924.39 

1.72028 

0. 745605 

5.8 

7732. 84 

0. 30982 

1.035939 

5.6 

j 5683.60 

1.64577 

0.753932 

7.8 

2612.81 

0.30803 

1.043459 

4.7 

: 53  '*9,  53 

1.54927 

0. 764236 

9.1 

2467.3] 

0.30455 

1.051602 

6.3 

: 5410.85 

1. 56000 

0.763330 

8.5 

2477.33 

0.30559 

1.051600 

4.6 

: 5104.37 

1.46075 

0.774234 

11.9 

2324.  93 

0. 30195 

1.060370 

5.4 

i 4810.89 

1.37208 

0.784642 

12.5 

2179.43 

0.29950 

1.069803 

4.2 

4544.00 

1.29984 

0. 793891 

15.9 

2047.01 

0. 29895 

1.079907 

5.4 

4249.49 

1.21429 

0. 805041 

19.3 

1901. 13 

0.29686 

1.090750 

4.7 

Temp. , 
degrees 


Vedam 

and 

Holton 


Pq  bars 


TABLE  III 


Fir  of  Hindi rz  and  Usual  Tail  Equations  to  Vedam  and  Holton  Data. 


Tbmlirz  Equation 

. _ v ....  i 


Usual  Tait  Equation 


\x  10  J 
bars  cc/gm 


Vffl  ee/gm  * x 10 
cc/gm 


B bars  C bars  cc/gm  Vjcc/gm  <rx  10 

cc/gm 


t 30 

8024 

2.7583 

0. 6600 

334 

2730 

0.2927 

1.0046 

178.0 

| * 40 

8934 

3.2767 

0.6398 

597 

2824 

0.2995 

1.0082 

161.0 

t 50 

8956 

3.2931 

0.6431 

599 

2833 

0.3004 

1.0124 

156.0 

| 60 

9004 

3.3443 

0. 6443 

618 

2853 

0.3037 

1.0173 

137.0 

g . 70 

8936 

3.3353 

0.6480 

657 

2822 

0.3046 

1.0228 

109.0 

| 80 

8760 

3.2747 

0.6535 

718 

2742 

0.3036 

1.0290 

81.5 
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Fig.  1 Temperature  dependence  for  pure  water  of 
the  parameters  PQ,  \ and  V»  for  the  TUm- 
lirz  equation. 

From  Table  I the  goodness  of  fit,  «■  • for 
the  Amagat  data  shows  a rather  good  uniform  fit. 
The  standard  error  is  about  the  same  as  the  quality 
of  the  data.  The  two  equations  fit  about  equally  well 
ao  that  a superiority  of  one  equation  over  the  other 
is  not  demonstrated.  These  flu  do  not  seem  to  dis- 
play any  inherent  limitations  in  the  applicability  of 
either  the  Tumlirz  or  UTE  to  pure  water  data. 


Fig.  2 Temperature  dependence  for  pure  water  of 
the  parameters  B,  C and  Vj  for  the  Usual 
Tail  Equation. 

The  very  accurate  data  of  Kell  and  Whalley 
is  more  revealing.  Excellent  fits  are  to  the  iso- 
therms obtained  below  90°  with  the  Tumlirz  equa- 
tion and  above  40°  with  the  UTE.  In  either  event, 
the  standard  error  is  about  5 ppm.  Figure  3 shows 
a plot  of  die  standard  error  at  each  temperature  for 
the  fit  of  both  the  proposed  equations  to  Kell  and 
Whalley  data. 


<r  IRPM) 


l 
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Fig.  3 Comparison  of  the  standard  deviations  <r 
of  fits  of  the  Tumlirz  equation  (left)  and  the 
Usual  Tait  Equation  (right)  to  the  Kell  and 
Whalley  d.ita. 


TABLE  IV 

Fit  of  Tumlirz  Equation  to  Pure  Water  Data  at  10°. 


Data 

P0  (bars) 

\bar 

cc/gm 

V*  cc/gm 

it  x 106  cc/gm 

Amagat  (1893) ! 

6228 

1858 

0.7019 

15 

Amagat  (1893)  II 

6570 

2038 

0.6900 

82 

Newton  & Kennedy 

6279 

1902 

0. 69-,  4 

22 

Kell  t Whalley 

6323 

1920 

0.6967 

6 

Vedam  & Holton 

6271 

1882 

0.7002 

26 

(low  pressure) 
Wilson  & Bradley 

6215 

1853 

0.7022 

24 

-5- 
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As  fits  were  made  to  data  from  different  in- 
vestigators, it  became  apparent  that  there  were 
large  differences  in  optimum  parameter  values  at 
the  same  temperature  as  illustrated  in  Table  IV  with 
the  TUmlirz  equation  for  pure  water  at  10°.  Exam- 
ination of  the  Tumlirzequation  shows  that  even  small 
changes  in  any  one  parameter  would  greatly  affect 
the  predicted  volumes  and  hence  the  quality  of  fit. 
However,  as  Tait  noted,  large  changes  in  one  par- 
ameter may  be  made  provided  compensating  changes 
are  made  in  the  other  parameters;  for  example,  in- 
creasing X and  PQ  while  decreasing  V*  would  have 
little  effect  upon  the  resulting  curve.  This  is  illus- 
trated in  Figure  4a  ( \ vs,  P0)  and  Figure  4b  (P0 
vs.  Vm  ) where  the  sloping  line  is  derived  from  an- 
alysis of  Kell  and  Whalley  data.  The  derivation  of 
these  compensating  relationships  is  explained  in  Ap- 
pendix II. 

In  comparingthefitsto  the  Kell  and  Whalley 
data  with  equations  (2),  UTE,  and  (3),  TUmlirz,  .wc 
see  that  even  though  DTE  parameter  B (Figure  2a) 
displays  a rather  smooth  behavior,  C (Figure  2b) 
shows  an  erratic  temperature  dependence  at  the 
higher  temperatures.  This  contrasts  with  the  smooth 
behavior  of  die  TUmlirz  parameters  (Figure  1)  over 
the  whole  temperature  region  except  at  50°.  Upon 
re-examining  the  <r  vs.  T dependence  for  the  two 
equations  in  Figure  3,  we  should  also  note  that  the 
precision  of  fits  of  the  Tumlirz  equation  vary  in  ac- 
cord with  the  errors  in  accuracy  estimated  by  Kell 


and  Whalley  as  a function  of  temperature.  On  this 
basis  we  select  the  Tumlirz  equation  to  represent 
the  PVT  properties  of  pure  water.  In  our  final  cal- 
culations we  have  included  Kell's**/  data  for  pure 
water  at  atmospheric  pressure 

Hie  TUmlirz  equation  and  the  temperature 
dependence  of  its  parameters  are  shown  in  Table  V. 
Comparisons  of  fits  obtained  with  our  equation  and 
that  observed  by  Kell  and  Whalley  are  shewn  in 
Table  VI.  It  is  quite  reassuring  to  note  that  the 
density  maximum  at  atmospheric  pressure  accord- 
ing to  our  results  occurs  at  4. 00°. 

Although  McDonald  does  not  consider  the 
Tumlirz  equation  as  a best  choice  for  representing 
the  properties  of  pure  water  we  believe  it  provides 
a good  base  for  an  analytic  function  to  represent 
specific  volumes  on  an  absolute  basis.  In  their  rep- 
resentations for  thePV  isotherm  data  both  McDonald 
and  Kell  and  Whalley  normalize  their  equations  to 
the  atmospheric  pressure  specific  volumes.  This, 
of  course,  makes  it  easier  to  include  later  and  more 
accurate  atmospheric  pressure  data.  It  also  makes 
calculation  of  thermodynamic  properties,  particu- 
larly the  thermal  expansion  coefficient,  more  in- 
volved. 


The  TUmlirz  equation  for  pure  water  also 
provides  a good  basis  for  representing  the  proper- 
ties of  sea  water  as  we  shall  see  next. 


i 


Fig.  4 Plot  of  X and  V*  vs.  P0  values  obtained  at 
optimum  fits  of  Tumlirz  equation  to  pure 
water  at  about  10  from  various  sources. 
The  straight  Une  is  the  compensating  rela- 
tionshipderived  from  Kell  and  Whalley  data. 
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Table  V 

(a)  Pure  Water,  0 - 100°,  8 ppm  and  100  - 150°,  15  ppm  fits  to 

Kell  and  Whalley  data  < X in  bars  cc/gm)  P,  P0  in  bars  and  14a  , V in  cc/gm. 

X = 1,788,316  + 21.55053  T - 0.4695911  T2  +3.096363  x lO^T3  - .7341182  x 10'5  T4 

P =*  5918. 499  + 58. 05267  T - 1.  1253317  T2  + 6. 6123869  x 10'3  T3  - 1. 4661625  x 10'5  T4 
o 

V»  = .6980547  - .7435626  10*3  T+  .3704258  x 10‘4  T2  -.6315724  x 10'6  T + . 9829576  x lO*8 
T4  - . 1197269  x 10'9  T5  + . 1005461  x 10' 11  T6  - .5437898  x 10‘14  T7  + . 169946  x 10* 16 
T8  - .2295063  x 10'19T9 


B.  Sea  Water 

The  most  extensive  sea  water  measure- 
ments have  been  made  bv  Wilson  and  Bradley  4// 
and  Newton  and  Kennedy8'. 

The  Wilson  and  Bradley  data  Is  published  in 
two  forms:  the  original  experimental  data  and  the 
smoothed  data;  the  isotherm  precision  to  be  + 20  ppm 
In  density.  With  a 22  term  parameter  equation  they 
obtained  a fit  with  a standard  deviation  of  80  ppm  and 
with  the  10  parameter  Tumlirz  equation,  which  they 
finally  chose,  they  obtained  a standard  deviation  of 
130  ppm.  Because  a Tumlirz  equation  was  used  as 
a part  of  the  smoothing  procedure,  the  smoothed 
data  cannot  be  used  to  test  the  applicability  of  the 


Table  VI 


Comparison  of  Standard  Deviations  of  Fits  to  Kell 
and  Whalley  Data  of  the  Equations  of  State  Kell  and 
Whalley  used  and  the  one  we  use. 


Temp 

Standard  Deviation 
Fisher  & Dial 

x 108  (cc/gm) 
Kell  & Whalley 

0.000 

9.6 

4,4 

10. 002 

8.0 

10.0 

19.  997 

7.4 

7.9 

24. 998 

7.0 

4.9 

24.  998 

5.8 

4.6 

25.005 

5.5 

4.3 

30. 1 12 

6.9 

11. 1 

39.  999 

8.3 

13.8 

50. 007 

6.7 

6.0 

60.001 

5.6 

4.3 

70. 003 

11.2 

13.6 

80. 003 

8.4 

4.2 

90.007 

6.7 

4.6 

100. 005 

10.2 

5.8 

110.005 

11.8 

10.3 

120.  007 

13,2 

9.6 

130. 010 

14.3 

5.8 

140.009 

18.5 

9. 1 

150.016 

22.8 

9.0 

Tumlirz  equation.  The  original  data  occurs  at  ir- 
regular temperature  intervals,  complicating  com- 
parison at  different  salinities.  Consequently  we 
have  used  polynomial  interpolation  to  provide  data 
at  even  5°  intervals;  it  is  this  interpolated  data  that 
will  be  referred  to  as  the  Wilson  and  Bradley  data. 

The  UTE  and  Tumlirz  equation  were  used 
to  fit  the  seawater  Isotherm  data  at  each  salinity. 
The  fits  were  on  the  order  of  30  ppm  for  both  the 
UTE  and  Tumlirz  equation;  therefore,  the  Tumlirz 
equation  was  selected  for  further  study  because  the 
quality  of  fit  does  not  restrict  the  selection  and  pro- 
vides a better  fit  to  pure  water. 

Figure  5 shows  x , PD  and  V®  vs.  salinity 
at  10°from  the  fits  of  the  Thmlirz  equation  to  Wilson 
and  Bradley,  and  Newton  and  Kennedy  data.  The 
parameters  do  not  appear  to  vary  in  any  consistent 
way  with  salinity  because  random  parameter  varia- 
tions resultingfrom  imperfect  data  obscure  the  sys  - 
tematic  variations  with  salinity. 


The  sea  water  data  was  normalized  to  pure 
water  by  constraining  X to  be  equal  to  the  value  for 
pure  water  at  the  same  temperature.  The  X values 
for  Kell  and  Whalley  data  were  polynomial  interpo- 
lated to  the  temperatures  of  the  sea  water  data  and 
least-square  fits  were  made  for  V®  and  P0  at  each 
temperature  and  salinity.  Figure  6 shows  the  re- 
sults at  10°  from  which  a linear  salinity  dependence 
is  now  evident.  A linear  variation  of  P0  with  salin- 
ity is  in  accordance  with  Tsunmann's  hypothesis 
which  states  that  Internal  pressure,  the  usual  inter- 
pretation for  P0,  should  increase  linearly  with  sal- 
inity. Also,  V®  is  seen  to  vary  linearly  with  sa- 
linity, Using  the  linear  variations  suggested  by 
Figure  6 we  may  writell/: 


Vp  =V®  - K^  + p +k  p 


(5) 


Equation(5)was  fit  to  the  Wilson  and  Bradley 
data  by  setting  V®  , X and  P0  equal  to  the  pure 
water  values  calculated  for  the  temperatures  at 
which  the  data  was  taken;  attempts  to  determine  the 
temperature  dependence  of  Kj  and  K2  were  not 
entirely  satisfactory  partly  because  the  data  were 
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obtained  at  varying  temperatures.  In  addition  it  v/as 
found  that  thermal  expansion  coefficients  calculated 
from  our  parameters  for  equation  5 did  not  agree 
with  the  recent  work  of  Bradshaw  and  Schleicher^/ • 
While  we  are  able  to  obtain  with  uur  oqua- 
tion  fits  to  the  Wilson  and  Bradley  data  with  stan- 
dard deviations  of  about  10"4  cc/gm  which  is  com- 
parable to  the  Fit  Wilson  and  Bradley  obtained  with 
their  equation  (1. 3 x 10*4)  we  feel  that  is  possible 
to  obtain  a better  equation  of  state  by  using  sound 
velocity  data  to  determine  I^and  thermal  expansion 
data  to  determine  Kj.  This  has  been  done  in  another 
paper15/. 


Fig.  5 Plots  of  PDand  V®  obtained  from  optimum 
fits  of  Tumllrz  equation  to  sea  water  data  at 
about  10°. 


DISCUSSION 
A.  Pure  Water 


Our  calculations  have  led  us  to  use  the 
Tumllrz  equation  to  represent  the  properties  of  pure 
water.  It  is  rather  Interesting  to  note  that  the  value 
of  V „ is  very  close  to  the  specific  volume  calcula- 
ted for  the  close  packed  form  of  water  In  most 
theories  which  treat  water  In  terms  of  a two  state 
model 1 . Along  another  line  of  thought  we  note  that 
David  and  Litovltz16/  calculate  the  temperature  de- 
pendence of  the  relaxing  and  non-relaxing  compres- 
sibilities for  water  and  find  a cross  over  point  in  the 
temperature  region  near  the  compressibility  mini- 
mum at  45  and  sound  velocity  maximum  at  72°.  In 
this  paper  we  have  seen  that  a crossover  in  the  pre- 
cision of  fits  of  the  Tumiirz  and  UTE  equations  also 
occurs  near  the  same  temperature. 

Since  the  open  packed  structure  (ice -like) 
predominates  at  lower  temperatures,  it  is  tempting 
to  associate  its  properties  with  the  Tumllrz  equation 
and  the  close  packed  structure  with  the  UTE.  We 
have  seen  that  the  UTE  equation  provides  a more 
precise  fit  to  Kell  and  Whalley's  data  at  higher  tem- 
peratures and  we  also  note  that  the  very  high  pres- 
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produce  a breakdown  of  the  open  structure  of  water. 
In  this  sense,  the  transition  from  the  Tumlirz  to  the 
UTE  equation  appears  to  be  consistent.  Using  the 
results  of  Davis  and  Litovltz  for  the  fractions  of 
open  and  closed  structure  we  attempted  to  fit  the 
Kell  and  Whalley  data  with  a linear  combination  of 
the  TUmllrz  equation  and  the  UTE.  The  results  so 
far  are  inconclusive  and  v.o  hope  to  pursue  the  mat- 
ter later,  based  in  part  on  the  work  of  Frank  and 
Quist.!!/ 

The  parameter  B in  the  UTE  and  P0  in  the 
Tumlirz  equation  have  been  referred  to  as  an  inter- 
nal pressure.  Since  no  experimental  data  as  yet 
exists  on  the  ultimate  tensile  strength  of  water,  we 
cannot  ascribe  an  interm olecular  significance  to  the 
internal  pressure.!®/  Going  back  to  the  two-state 
model  we  note  that  the  lower  internal  pressure 
values  of  B would  be  attributed  to  the  disordered 
close  packed  state  and  the  much  higher  P0  values  to 
open-packed  state.  Again  this  seems  to  be  reason- 
able from  a qualitative  viewpoint. 


I 

i 


f 

I 

I 

i 

& 

i 


Slfmity  S (%•) 


Fig.  6 Plots  of  P0  and  Vm  obtained  from  optimum 
fits  of  Tumlirz  equation  to  sea  water  data  of 
Wilson  and  Bradley  at  10  when  X is  con- 
strained to  be  the  value  for  pure  water. 


sure  data  of  Vedam  and  Holton10/  at  lower  temper- 
atures shows  a lower  standard  deviation  when  the 
UTE  is  used  for  fitting. 

Although  the  accuracy  of  the  data  at  higher 
temperatures  and  very  high  pressures  is  notasgood 
as  that  of  the  low  temperature  Kell  and  Whalley 
data,  the  fact  that  the  UTE  provides  better  a fit  at 
the  high  temperature  and  pressures  may  be  signifi- 
cant. Both  increasing  temperature  and  pressure 


B.  Sea  Water 

Although  we  have  not  established  the  tem- 
perature dependence  of  the  parameters  Kj  and  Kj  to 
our  satisfaction  in  this  paper,  we  have  clearly  dem- 
onstrated the  linear  salinity  dependence  of  the  cor- 
rections to  Vco  and  P0  as  shown  in  equation  5.  This 
promises  to  be  a very  useful  equation  in  that  atmo- 
spheric pressure  data  on  the  density,  sound  ve- 
locity, and  heat  capacity  can  be  used  to  determine 
the  high  pressure  properties  of  electrolyte  solu- 
tions. Furthermore,  it  promises  to  be  useliil  for 
determining  partial  molal  properties  of  solutions  at 
infinite  dilution  from  data  at  high  concentrations. 
Unfortunately,  most  data  at  atmospheric  pressure 
and  elevated  pressures  up  to  now  on  solution  densi- 
ties are  not  accurate  or  precise  enough  to  do  a de- 
finite study  utilizing  equation  5 in  this  matter.  12/ 
Experiments  are  in  progress  at  various  labora- 
tories which  promise  to  yield  data  of  the  necessary 
accuracy.  It  may  also  be  possible  to  calculate  the 
properties  of  solutions  with  more  than  one  solute  by 
additions  of  their  respective  values  of  Kj  and  K2. 

The  principal  difference  In  the  form  of  equa- 
tion 5 from  that  used  by  Eckart  and  Wilson  and 
Bradley  is  that  X is  not  a function  of  salinity.  This 
may  be  an  artificial  constraint.  More  pre- 
cise solution  data  will  determine  the  utility  of  this 
equation.  Yayanos  11/  has  used  another  equation 

PV4  = A - V/V  (6) 

o 

to  analyze  his  PV  isotherms  for  water  and  NaCl.  In 
further  work  we  expect  to  test  this  equation  also. 
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intervals  about  the  optimum  value  found  by  the  pre- 
vious search.  This  search  was  coded  intoacomputer 
routine  which  terminated  when  the  optimum  B was 
localized  to  + . 001  bar. 

The  fitting  procedure  for  the  Tumlirz  equa- 
tion is  quite  similar.  Values  of  VQ,  k , and  cr  are 
calculated  for  selected  values  of  PQ.  A search  for 
the  optimum  is  continued  until  PQ  is  locallzed+  ,001 
bar. 

By  computerizing  nonlinear  least-square 
calculations,  we  are  able  to  obtain  the  optimum  fits. 
In  all  cases  the  standard  error  is  computed  to  pro- 
vide a simple,  uniform  measure  of  the  quality  of  fit; 
indeed,  the  minimization  of  this  quantity  is  the  es- 
sence of  least- square  fitting. 


APPENDIX  I: 

COMPUTATIONAL  PROCEDURE  FOR  NONLINEAR, 

LEAST-SQUARE  CURVE  FITTING 

Nonlinear  least-square  curve  fitting  re- 
quires the  use  of  iterative  approximation  to  calculate 
parameters  of  the  curve.  We  use  the  fit  of  the  UTE 
equation  (2),  to  the  Kell  and  Whalley  data  at  10.002® 
to  illustrate  the  iteration. 

We  define  <r , the  standard  error,  by  equa- 
tion (4).  Least  square  curve  fitting  implies  the  se- 
lection of  parameters  Vj,  C,  and  B to  minimize  a , 
Of  these,  B is  nonlinear. 

If  the  value  of  B in  equation  (2)  were  known, 
we  could  use  conventional,  linear  least-square  fit- 
ting to  determine  V^  and  C and  calculate  the  re- 
sulting <r . Table  VII  shows  the  results  of  these  cal- 
culations for  tabulated  values  of  B.  It  is  readily 
seen  that  <r  smoothly  comes  to  a minimum  value 
near  B = 2900  bars.  The  nonlinear,  least-squane  fit 
is  implemented  by  searching  successively  smaller 


APPENDIX  11: 

CORRELATED  PARAMETER  ESTIMATES 

Parameter  estimates  obtained  from  dif- 
ferent sets  of  data  show  both  random  and  systematic 
variations.  Experience  with  the  Tumlirz  equation 
has  shown  that  there  are  large  variations  in  param  - 
eterestimates  obtainedfrom  all  but  the  most  pre- 
cise data.  Figure  5 shows  that  random  variations  of 
parameters  for  seawater  data  are  so  large  that  they 
obscure  any  systematic  variations  with  salinity. 

By  considering  different  sets  of  data,  we 
have  found  that  the  estimates  of  the  parameters  axe 
correlated  as  shown  In  Figure  4.  Also,  any  one 
parameter  can  be  changed  without  significantly  de- 
grading the  quality  of  fit  provided  that  compensating 
(i.e. , the  correlated)  changes  are  made  in  the  other 
parameters.  Table  V maybe  considered  to  tabulate 
the  values  of  X and  \m  necessary  to  compensate 
for  variations  of  P0.  Indeed,  when  the  X vs.  P0  and 
V»  vs.  P0  columns  from  Table  V are  plotted  in 
Figure  4 the  previously  obtained  parameterestimates 
are  seen  to  agree  closely. 


Table  VII 


Least-Square  X , V0,  and  <r  from  Fit  of  Tumlirz 
Equation  to  Kell  and  Whalley  Data  at  10°C. 


Pn(bar) 

lambda(bar-cc/gm)  Vn(cc/gm) 

<r  (ppm) 

5600 

1,530.894 

0.727055 

60.4 

5800 

1,634.105 

0.718665 

42.5 

6000 

1,740.662 

0.710276 

26.0 

6200 

1,850.564 

0.701888 

11.2 

6323. 46 

* 1,920.082 

0.  696711 

6.2 

6400 

1,963.813 

0. 673502 

8.4 

6600 

2,080.406 

0.685117 

20.7 

6800 

2.200.346 

0. 676733 

33.7 

’optimum  value  found  after  continued  search 


Similarly  correlated  parameter  estimates 
occur  with  the  UTE,  These  correlations  explain  the 
success  of  previous  work  which  constrained  the  C 
parameter  because  B values  could  then  be  selected 
Which  largely  compensated  for  the  non-optimum  C 
values  chosen21^. 

The  existence  ofthe  compensating  relations 
also  provides  a rationalization  for  constraining  a 
parameter  of  the  seawater  equation.  We  know  such 
constraints  won't  degrade  the  fit  and  that  the  com- 
pensating relations  will  reduce  the  erratic  variations 
of  the  remaining  parameters.  In  particular,  con- 
straining the  X parameter  to  its  value  for  pure 
water  at  the  same  temperature  produced  the  linear 
variations  of  V*  and  P0  with  salinity  which  Tamman 
had  suggested. 
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